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1 A two dimensional ruin problem 

In this paper we consider a particular two dimensional risk model in which 
two companies split the amount they pay out of each claim in positive pro- 
portions 6i and 62 with 61 + 62 = 1, and the premiums according to rates ci 
and C2- Thus, the risk process f/j of the z'th company satisfies 

Ui{t):=-6,S{t) + Cit + Ui, 2 = 1,2, 

where Ui are the initial reserves. We will work with a spectrally positive Levy 
process S{t), that is Levy process with only upward jumps that represents 
the cumulative amount of claims up to time t. In particular we focus on the 
classical Cramer-Lundberg model: 

Nit) 

sit) = j2^k, (1) 

fc=l 

where N{t) is a Poisson process with intensity A and the claims Ck are i.i.d. 
random variables independent of N{t), with distribution function F{x) and 
mean E[ak] = fi~^- We shall assume that the second company, to be called 
the reinsurer, gets smaller profits per amount paid, i.e.: 

di 62 

As usual in risk theory, we assume that pi > p := -, which implies that in 
the absence of ruin, Ui{t) ^ 00 as t —>■ 00 {i = 1,2). Ruin happens at the 
time r = t{ui, U2) when at least one insurance company is ruined: 

r(Mi,M2) :=inf{t> : f/i(t) < or f/2(t) < 0}, (3) 

i.e. at the first exit time of (?7i(t), U2{t)) from the positive quadrant. In this 
paper we will analyze the perpetual or ultimate ruin probability: 

tp{ui,U2) = P[t{ui,U2) <oo]. (4) 

Although ruin theory under mult i- dimensional models rarely admits analyti- 
cal solutions, we are able to obtain in our problem a closed form solution for 
(jlj) if ai are exponentially distributed with intensity /i. 




> Ml 



Figure 1: Geometrical considerations 



Geometrical considerations. The solution of the two-dimensional 
ruin problem (jl]) strongly depends on the relative sizes of the proportions 
6 = (61,62) and premium rates c = (ci,C2) - see Figured! If, as assumed 
throughout, the angle of the vector 6 with the ui axis is larger than that of 
c, i.e. 62C1 > 61C2 we note that starting with initial capital (mi, U2) G C in the 
cone C = {(^1,^2) : U2 < {62/Si)ui} situated below the line M2 = ((^2/5i)mi, 
the process (f/i, U2) ends up hitting at time r the Ui axis. Thus, in the domain 
C ruin occurs iff there is ruin in the one-dimensional problem corresponding 
to the risk process U2. 

One dimensional reduction. A key observation is that r in ([3]) is also 
equal to 

t{ui,U2) = inf{t > : S{t) > b{t)}, 

where b{t) = m.m{{ui + Cit)/6i, {u2 + C2t) / 62} ■ The two dimensional problem 
dlj) may thus be also viewed as a one dimensional crossing problem over a 
piecewise linear barrier. 

In the case that the initial reserves Ui and U2 are such that {ui,U2) G C, 
that is, U2/S2 < Ui/61, the barrier b is linear, b{t) = {u2 + C2t)/(52, the ruin 
happens always for the second company. Thus, as we already observed, the 
problem (jlj) reduces in fact to the classical one-dimensional ultimate ruin 
problem with premium C2 and claims ^20", i.e. 



^{ui,U2) = ^'2(^2) := ^(7-2(^2) < 00) 



where T2{u2) = inf{t > : U2{t) < 0} and 4'2{u2) is the ruin probabihty of 
U2, with 1/2(0) = U2- For the model ([T]) the Pollaczeck-Khinchine formula, 
well known from the theory of one-dimensional ruin (see e.g. [8] or [1]), 
yields then an explicit series solution for ?/;(Mi,ti2) = ^2(^2) in the case of 
a general claims distribution. For the phase- type claims {(3,B), i.e. with 
P[cr > x] = fBe^^l, the ruin probability may be written in a simpler matrix 
exponential form: 

^2(^2) = rie^2HB+br,)u,-^ 

with 77 = —(3{—B)^^ (see for example (4) in f2j), and in the case of expo- 
nential claim sizes with intensity /i, it reduces to: 

^2(^2) = C2e-(^^/'^)"^ , (5) 

where 72 = ^u - X52/C2 = /i - \/p2 and C2 = g = ^. 

The rest of the paper is devoted to the analysis of the opposite case, 
^2/^2 > ui/^i and is organised as follows. Section [2] is devoted to the 
Laplace transform of the ruin time in the case S is of the form ([T]) with 
exponential jumps. Subsequently, in Sections [3] and H] we derive the Laplace 
double transform in (mi, M2) of the ruin probability ^(mi, M2) US is a general 
spectrally positive Levy process. Finally, in Section [5] this Laplace transform 
is explicitly inverted, in the case of exponential claim sizes. 



2 The differential system for the exponential 
claim sizes case 

In this section we provide a system of partial differential equations for the 
Laplace transform 

^Piuu U2, s) := E[e-^^("i'"^)l{,(,„„,)<oo}] (6) 

of the ruin time t{ui,U2) in ([2]) in the case that S is given by a com- 
pound Poisson process ([T]) with intensity A and with claims sizes ai that 
are exponentially distributed with parameter /x. The memoryless property 
of interarrival times and claim sizes opens up the possibility of embedding 
the (discontinuous) Cramer-Lundberg processes (f/i,f/2) into a continuous 
Markov-modulated fluid model. Informally, this is achieved by a transfor- 
mation that replaces the jumps of (f/i,f/2) by a linear movement in the 



direction {—61,-62) of duration equal to the size of the jump, creating 
thereby a new continuous semi-Markovian model (f/i,f/2) called the fluid 
embedding of {Ui,U2)- As the process (f/i,f/2) crosses boundaries continu- 
ously and has exactly the same maxima and minima as (f/i, f/2), first pas- 
sage problems may be easier to handle for {Ui,U2) than for the original 
process {Ui,U2). A formal construction of the fiuid-embedding is given in 
the Appendix. If we write f{ui,U2) = inf{t > : m.m{Ui(t) , U2{t)} < 0} 
for the joint ruin time of {1/1,1/2), it follows from the definition of {1/1,1/2) 
that t{ui,U2) = I{f{ui,U2)), where I{t) denotes the time up to time t that 
{1/1,1/2) was increasing. In particular, 

^{ui,U2,s) = ^[e-^'("("^'«^»l{,(„„„,)<oo}|ro = 1]. (7) 

Setting (t){ui,U2,s) = ^[e~*^(^("i'"2))l|-(„^ ^^^^^^jIFq = -1] and writing ^„, 
and (pui for the partial derivative of ip and with respect to Ui we have the 
following characterization of and ip. 

Theorem 1 For 62U1 < 61U2 it holds that {il){ui,U2,s), (f){ui,U2,s))~'' solves 
the Feynman-Kac system: 

ci \ f ipm\ f C2 \ f ipu2 
-6i J\(l)u, J^\0 -62 J \ <pu, 

-X-s X \ f ip \ f 



' A* ~t^ ) \ ^ ) \^ 

with the boundary condition: 

i^{ui,fui) =C2e-''^*"^ forallu^y^, ^g^ 

0(0, U2) = 1 for all U2 > 0. 

Proof. Conditioning on the first shared claim occurrence epoch we obtain: 

■0(-ui, U2, s) = e~^^e~''^il){ui + cih, U2 + C2h, s) 

+ / Ae~^* dt / /ie~^''e"''V('"i + Cit - 6iz, U2 + C2t - 62Z, s)dz 

Jo Jo 

nh POD 

+ / Ae"^* dt / fie-^'dz; (9) 

Jo J^<i+£it^«2+£2i 

"1 *2 

(j){ui,U2, s) = / fie~^^^lj{ui — 6iZ,U2 — 62Z, s)dz (10) 

Jo 

5 



Both integrals on the LHS of ([H]) go to as /i — > 0. This imphes that ^^ is a 
continuous function with respect to Ui and U2- Note that Hi±£i* /\ m±£2t _ 
"^j"^^* since we hve in the upper cone. This gives after simple manipulation: 

ip{ui + Cih,U2 + C2h,s) -ip{ui,U2,s) e-^^e-''^ -I , u \ 

; \ ; W[Ui + Ci/l, U2 + C2h, s) 

n n 

1 rh i'{ui+cit)/5i 

+- / Ae"'^* dt / /xe"'^^e"''V('"i + Cit - 6iz, U2 + C2t - 52Z, s)dz 



h Jo Jo 

-h 



hJo 



Note that the last 3 terms on the LHS of above equation have limits as 
/t — ^ 0, since ?/' is a continuous function. Thus if) is a, different iable function 
of {ui,U2). Taking /i — > we derive 

Ci^„^ (Mi, ti2, S) + C2^U2 (Mi, W2, s) + (-A - s)V'(Mi, ti2, s) 

+A / ' ' /ie-^^V(Mi - (^i-z, M2 - 52Z, s)dz + Ae-^"^/^^ = 0. 
Jo 

In view of ( TTU]) this equation is equivalent to 

CiV'm («i> "2, s) + CzV'^j (mi, M2, s) + (-A - s)V^(mi, M2, s) + A0(Mi, M2, s) = 0. 

This gives the first equation in the Feynman-Kac system. To derive the 
second one, we apply integration- by-parts formula to (ITOl) : 

f^ d 

0(mi,M2,-s) = / --{-e'^'')ij{ui-SiZ,U2- S2Z,s)dz 
Jo dz 

= i;iu,-S,z,U2-52Z,s)i-e-n\^ 

f°° _ d 
+ e ^''—iIj{ui — 6iZ,U2 — S2Z, s)dz 



dz 



00 



ip{ui,U2,s) - 61 / ipuiiui - 6iz,U2- 62Z,s)e ^""dz 





00 



- S2 'ipuiiui - 6iz,U2 - S2Z, s)e ^""dz 
Jo 

= t/j{ui,U2,s) - (5iyU"V«i(Ml,M2,s) - 62fJ.'^(f)u2{ui,U2,s), 

which gives the second equation in the Feynman-Kac formula. The boundary 
conditions follow immediately. D 
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The above system may also be reformulated as a second partial-differential 
equation in terms of ■ip{ui,U2, s) only. To that end, we define a linear trans- 
formation (x, of (V^5 0) by x{r, w, s) = ^(x, y; s) and ^(r, w, s) = 0(x, y; s) 
where {x,y) = {x{r,w),y{r,w)) are given by 

X \ ^ f 5i ci\ f r 
y J \52 C2 ) \w 

with inverse transformation: 

r \ _ ,-1 ( C2 -Ci \ ( X 



w J \ -S2 Si J \ y 

where d := 61C2 — S2C1 < 0. In the next result a PDE is derived for x- 
Corollary 1 The function 

h{r,w) ■.= €^'^6-'^^+'^''" x{r,w;s) 

solves the equation 

hrw + /iA /i = (11) 

with the boundary conditions: 



h{r, 0) = C2e~(^2-M) r j^^ ^11 ^ > Q^ 

/i^(r,-^) = -Ae^'^"^^+'^^ forallr>0. 



Proof: Note that: 

• during drift periods, r is constant and w increases, at rate w' = 1 

• during jump periods, w is constant and r decreases, at rate r' = —1. 

Thus, in these new coordinates, the time is split between moving into 
the direction of the axes and moving away from the axes. In particular, 
at any time T we have T = T^ + Tr, where T^j/Tr are the total times of 
growing reserves (upward drifting) /shrinking reserves (jumping). Note that 
w = Wo + Tw and r = tq — % hold. 

In terms of the (r, w) variables, the Feynman-Kac system becomes: 



with 

X{r,0,s) =C2e-^^^ 

e(r,-|r,.) =1, ^''^) 

for all r > 0. Recall that the upper cone is described be inequalities r > 
and w <0. Following the steps used in the proof of differentiability of ip one 
can prove that x is in class C^. 
Eliminating ^, we find 

Xrw - (A + s)xr + jJ^Xw - snx = (14) 

with 

x(r,0,s) = C2e-^2^ 

\-'{{X + s)x{r, -|r, s) - x^r, -|r)} = 1. ^'^^ 

We may remove the linear terms by switching to the function h in terms of 
which we get the stated result. D 

3 Probabilistic solution 

One way to obtain the Laplace transform of the joint ruin probability is to 
solve the above systems numerically. Here, we pursue a different approach, 
by establishing first a general analytical representation of the solution that 
holds for a general spectrally positive Levy process S. 

Noting that the process {Xi,X2) = (t^i/^i, t/2/52) has the same ruin 
probability as the original two-dimensional process {1/1,1/2), we can restrict 
ourselves without loss of generality to the process (Xi, X2). In the sequel we 
will write ip{xi,X2) for the joint ruin probability ([3]) - (jl]) corresponding to 
the process (Xi,X2). 

Proposition 1 If X2 > Xi, then 

POO ^ 

^[xx,X2):=\-^[xx,X2)= \ V^2(^)-P(xi,T)(dz), (16) 

where 







T = T(xi,X2) = ^^^, (17) 

Pi -V2 



P(,.,,T)(d^) = p., ( infXi(s)>0,Xi(T)Gciz 
with Px denoting P conditioned on {Xi = x}. 



Proof: In view of the definition (j3]) of r we see that 

i^{xi, X2) = P (min{Xi(t), X2(t)} > for all t > 0) , 

where (Xi(0),X2(0)) = (xi,X2). Next, we note that, if X2 > xi, it holds that 
the minimum 

min{Xi(t),X2(t)} = minjxi - X2 + (pi - ^2)^, 0} + X2(t) 

is equal to Xi[t) for t < T and X2{t) for t > T, where T was defined in 
flT7|l . We have also Xi{T) = X2{T). Subsequent application of the Markov 
property of X2 at time T shows that 

^(xi, X2) = P(,.,,,,)(Xi(t) > for t < T, X2{t) > for t > T) 

= j P,, (x,{T) e dz, inf Xi(s) > o\ P, (inf ^2(5) > o\ 

= j P,, (x,{T) G dz,MX,{s) > o\ i^^iz). 

D 

In Section m we obtain the double Laplace transform of ip{xi,X2) in 3^1,3^2, 
which we invert in Section [5], in the case of exponentially distributed jumps, 
using Bromwich type contours. 

4 Double Laplace transform in space 

Let S{t) now be a general spectrally positive Levy process and denote by 
Ki{9) the Laplace exponent of the spectrally negative Levy process Xj(t) = 
Pit- S{t), 

E[e^^'W] = e"»(^)*, z = l,2. (18) 

We may obtain directly the double Laplace transform in space of ip{xi,X2), 
by exploiting for X2 > xi the integral representation in Proposition [1] and 
for X2 < Xi the explicit formula of the Laplace transform in x of the one- 
dimensional ultimate ruin probability 'ip2{x)- We will use the following re- 
sults: 



1. If Kj(0+) > 0, the Laplace transform with respect to the starting point 
of the ultimate survival probability ipi{x) = Px{mit>Q Xi{t) > 0) (see 
e.g. Bertoin ^, Thm. VII.8): 

i^.ne) := / e-^>,(x) dx = ^^ z = 1, 2. (19) 

Jo i^i{(^) 

2. The resolvent of a spectrally negative Levy process killed as it enters 
the nonpositive half-line, due to Suprun (see also Bertoin [5], Lem. 

1]): 

f 6""*^,, (inf Xi(s) > 0,Xi(t) G dz\ dt 

= [exp{-q+{q)z}W^'^\x,) - ^,,>,}W^'^\x^ - z)] dz, (20) 

where q^{q) largest root of Ki{a) = q and W^''^ : [0, oo) -^ [0, oo) is a 
continuous and increasing function (called the q- scale function oi Xi{t)) 
with the Laplace transform: 

e-»^W^i)(y)dy = Ma) - q)-\ a > g+(g). (21) 

Now we obtain the double Laplace transform of the non-ruin probability with 
respect to the initial reserves: 

/"OO /-oo 

Jo Jo 
Note that 

V'(P;'?)= / / e"^^^e"'^^^?/'(xi,X2)dx2dxi 
Jo Jo 

/•CXD /"OO 

+ / / e-P^ie-''^2^(a;i,X2)dx2dxi. 

Jo Jxi 

The first Laplace transform is given by 

oo PXl -| 

/ e-P'''e-''''^^Jx2)dx2dxi = -(V^2)*(P + q) ■= ^• 
Jo P 
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Writing s = p + q and r = (pi — ^2)0' we see from fll6p and fl20p that the 
second Laplace transform is given by 



00 ^00 



Jxi 



e ^^'e '^'''^i){xi,X2)dx2dxi 



'0 

Pi - P2 r,- 

~~r~\ — K^2 

Ki[s) — r 



-SXl 



dx, / ^2(z)[e-''^MwM(xi) - l|,<,,}iyM(xi - z)]dz 



nqHr))-i^,ns)]:=C-B, 



where for the calculation of quantity B we used ( 1191) and (1211) . In view of 
f fT9|) we note that the quantity A — B is equal to 



{^,ns)K2is) 



4(0+) 



(22) 



p(Ki(s)-r) p{Ki{s)-r)' 
Similarly, since K2{6) = i^i{d) + (P2 — Pi)0, we see that C can be written as 

4(o+)(pi-P2) 4(o+)(pi-P2) 

«2(g+(r))(fi:i(s) - r) [fi:i(g+(r)) + (p2 - Pi)g+(r)] (/«i(s) - r) 
4(0+)(pi-p2) 
(ki(s) - r)(r + {p2 - Pi)g+(r)) ' 

Putting everything together we find: 

Proposition 2 The double Laplace transform ip is given by 

[<(0+) + {P2 - Pi)][r + {pi ~ P2){p - q^{r))] 



ip{p,q) 



P[r + {P2 - Pi)q^{r)] {ki{s) - r) 



4(0+) 



p{Ki{p+q) -q{pi-p2)) 



1 + 



P 



g-g+(g(pi-p2)) 



(23) 



4.1 Exponential claims 



In this section we specialize the above result to the classical model ([T]) where 
the jumps are exponentially distributed with parameter /x (cij ~ E(Ai)) and 
we write pi = Ci/6i, i = 1,2. In this case the characteristic exponent of Xi is 
given by 

Xa 

Ki[a)=pia ■ — , 1 = 1,2. 

fi + a 
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In particular, in view of the form of ki and k,2 it can be verified that k,2{0+) 
P2 — p (with p = A//i) and k{p, q) = \og E[eP^'^^-^^^'^^^^^^] is equal to 

ifi+p + q) 
where r = {pi — p2)q and s = p + q and 



,. -{P2q + Pi{q + 7i)) - \/(P2g + Pi{q + 7i)) - ^PiqP2{q + 12)^.. 

.X -{P2q + Pi{q + 7i)) + V(P2g + Pi{q + liW - ^PiqP2{q + 72) 
"^(^) = 2^^ ' 

with 'ji = fi — X/pi, i = 1,2. For later reference we note that 
^i(O) = -71 
^i(-72) = ^(^-py,Z2i-l2) = ^(^-pY (25) 

P2 \Pl J P2 \Pl J 

with x~ = niin(x, 0), x~^ = max(x, 0). Noting that q~^{q{pi—p2)) is the largest 
root of Ki(a) = q{pi—p2) and -22(5') is the largest root of Ki{v + q) = q{pi—p2) 
we identify 

g+(g(pi -P2)) =Z2{q) + q. 

In view of fl23|) we thus arrive at: 

Corollary 2 If S is given by (QP with at ~ E(/i), then ip is given by 

7f X if^ + P + q)iP2 - p) .r>a\ 

PPi{zi{q) -P)z2{q) 

5 Spectral representation 

In this subsection we invert the Laplace transform fl2Bl) of the ruin probabil- 
ity ip for exponential claim sizes. To perform the inversion we shall employ 
the method of residues. For an overview of the theory of Laplace transforms 
and complex analysis see e.g. Widder [9] or Ahlfors [3J. The method of 
residues leads to an explicit analytical representation of the survival proba- 
bility 'il){xi,X2) given in the following theorem. 
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Theorem 2 Let X2 > Xi and let S is given by (|7]j with ctj ~ E(/i). Then it 
holds that 

T, , jl-C,e-''^^+uj{x^,X2), ^fp<i> 

W(Xi,X2) = < ^^ 

[ 1 - Cie-^i^'i - C2e-^2X2 ^ 22e-73Xi-72X2 + tu(xi, X2) e/se, 
where 73 = — ( p — — ) Q'^'i 
^■(^^^^j _ P2-P r- ^xiafa)+x2g /(g) sm(Kg)a^i) + Kq) cos{b{q)xi) ^^ ^^7) 

g± = (Va± VPlP)^ 

P1-P2 

/(g) = (p + g + a(g)) and 

, . -(PiP- A + pag + Pig) 

a(g) = 



biq) 



2pi 

A/4pi(p2gp + P29^ - Ag) - (pifi - A + P29 + Pig)^ 

2pi 



To prove this resuh, first observe that '?/'(xi, X2) can be recovered from ip{p, q) 
using Mellin's formula, as folows: 

— / / i>ip,q)e^'''e^'Pdpdq, (28) 

where a > 0. The next step consists in iteratively evaluating this double 
integral (first w.r.t. p and then w.r.t. q) using Cauchy's theorem. The result 
of the first inversion is given in the next result: 



Lemma 1 For a > and fixed q with 3?(g) > it holds that 

— / ^(p, q)e^^^dp = ^^ - e^^^^^'^^giq), (29) 

where g{q) is given by 



qiq) = '^P2~^)'^^ + ^i'^^) + ^) (30) 

q{fip2 - A + p2q) 
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Figure 2: Bromwich contour 



The proofs of results that are not developed in the text can be found in the 
Appendix. 

In view of equations (iT9l) and ([5]) we recognize the first term in ( !29l) as 
the Laplace transform of '?/'2(^) — ^ ~ C2e~'''^^- The inversion of the second 
term relies on the following properties of g{q) and Zi{q), that were defined in 
fl30|l and (^ respectively: 

Lemma 2 (i) The functions g and Zi, are analytic in the set 

g = {gGC:g^{0,-72}U[g+,g_]}, 

where —72 > q- and qg{q) remains bounded if \q\ —>■ 00. 

(a) Let q^ = q ± ie with q G [g+, g_] and e > 0. // e | 0, then 

zi{qt) -^ z'iQ) ■= 0(9) - iKl)^ ^i(C) -^ z+{q) := a{q) + ib{q). (31) 

In order to ensure that we can calculate the inversion of the second term in 
f l2U]) using the method of residues we fix < a < 72 and replace g{q) by 
9{'1)/{q + o,) (note that in view of Lemma [2] the latter is 0{q~^) as |g| — > cx)). 
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Denote by / and ha the Laplace inverses of g{q) and g{q)/{q + a), that is 

/•OO /'OO px 

Jo Jo Jo 

Then it follows that / can be recovered from h^ by 

f{x)=lim^h^{x). (32) 

alO dx 



To complete the inversion of ■?/'(p, q) we are thus led to evaluate the integral 

r"Q+ioo 

2Txi 



1 pa+ico r \ 

— / ka{q)dq where ka{q) := J^ilLg^'^s+^i^iW. (33) 

^7r« Ja-toc^ q + a 



Using Lemma [2] we choose a Bromwich contour F/j ,. that encloses the poles 
of e^^^'^^^^g{q) while the cut of the square root zi{q) is not enclosed (see 
Figure [2] and the proof of Lemma H] for formal definition of F^ ,.). We note 
that this is a standard approach to calculate integrals of the form fl33l) (e.g. 
Ahlfors [3] or see Pervozvansky [6] for a recent application to the calculation 
of one- dimensional ruin probabilities). 

Recalling from Lemma [2] that g{q) has two simple poles, in g = and 
q = —72, Cauchy's theorem implies that 

— - (h ka{q)dq = ReSq=oka{q) + ReSq=_^2/ca(g) + Resq=-aka{q)- (34) 

The next step consists in evaluating the residues in fl34l) . which is a matter 
of straightforward calulations: 

Lemma 3 Writing C2 = C2 + ^^i^^^ it holds for a > that 

Resg=.aka{q) = g{-a)e-^'^+^^'^^-''\ Res,=oka{q) = ^e'^^^^ 



Resq=_^,kaiq) = ^l^e^i (-^2)^1 --^2^2 _ 



a -72 
Next we turn to the left-hand side of the formula (E 
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Lemma 4 For a > it holds that 

lim lim — - (h ka{q)dq = -— / ka{q)dq - uJa{xi,X2), 

40 R^oo 2m Jy^^ Z-ni Ja-ioc 

where u!a{xi,X2) is given by [2^ but with dq replaced by dq/{q + a). 

Proof of Theoreml^ln view of (152]) . (1341) and Lemmas [3] and HI the final 
result is obtained by first differentiating the residues in Lemma [3] and Ua with 
respect to X2 and subsequently letting subsequently a tend to zero. Taking 
note of the facts 

d 

lim2;i(— a) = Zi(0) = —71, limag(a) = C2 and lim- — LJa(xi,X2) = uj(xi,X2) 

aiO aiO aiO 8X2 

completes the proof (where the latter follows using the dominated conver- 
gence thoerem). D 



A Appendix 

A.l Formal construction of the fluid-embedding 

A formal construction of the process (f/i, U2) is as follows. Let ti,T2, . . . and 
ai, (72, . . . denote the subsequent inter-arrival times and claim sizes. Note that 
these form sequences of i.i.d. exponential random variables with parameters 
A and /i respectively. Define the switch times Sn by 6*0 = and, for n > 1, 

S2n-1 = S2n-2 + ^nj S2n = S2n~l + (^71 if ^0 = 1 
S2n-1 = S2n-2 + <^ni 8211 = S2n-1 + Tn if ^0 = ~1 

and construct Y taking values in { — 1, +1} by setting 

Ys^ = -Ys^_, forn > 1 with Fq G {-1, +1}. 

Then F is a two-state Markov chain indicating whether {Ui, U2) is increasing 
(state -|-1) or decreasing (state —1); more precisely, denoting by 

Ht) = [ i{n=i}ds 

Jo 
the total time up to t that Y has spent in state +1, we set 

Uiit)=pilit)-S,it-Iit)), 2 = 1,2, 

and the construction is complete. 
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A. 2 Proof of Lemma [T] 

By performing a partial fraction decomposition (in p) it follows that 



^{p,q) 



P q{fJ'P2 - A + qp2) p - zi{q) q{fip2 - A + gpz) 
14(0) 1 



pK2{q) p-zi{q) 



9{q)- (35) 



Since f^ e '^^e'^^dt = {p — c) ^ for p > c, the result follows by inverting ( l35ll 
term by term. D 

A. 3 Proof of Lemma [2] 

(i) Noting that the argument of the square root in ( 12^ is positive if q is real 
and g < g+ or g > g_, it follows that zi{q) is analytic outside the cut [g+, g_]. 
Since, furthermore, the denominator of g has no roots in Q, we see that g{q) 
is analytic in the set Q. The asymptotics directly follow from the form of g. 
(ii) Employing the standard definition of the square root z ^^ \fz (with 
the cut along the negative half-line) and appealing to the definition of z\ 
and the continuity of the argument Arg(2;) and modulus \z\ imply that the 
convergence in (|3T1) holds true. D 

A. 4 Proof of Lemma [4] 

Consider the contour F/j ,. that is given in Figure [2l i.e. Fr ,. consists of the 
line segments \a — iR, a + iE] , [g_ + e — ze, g_ + e + ze] , [g+ — e — ie, g_ + e — ze] 
and [g+ — e + ie, g_ + e + ie\ and two quarter circles in the left half-plane 
joining g+ — e + ie and a + iR, and, g+ — e — ie and a — iR, respectively. 

By taking the limits of i? — >• oo and subsequently letting e | and using 
that the integrals of the quarter-circles tend to zero (in view of the fact that 
9{q)/{q + (^) = 0{q^'^) as |g| — > oo, cf. Lemma[2](i)) we find that the contour 
integral ^ ^p ka{q) dq converges to 

— - / ka{q)dq+-— ka{q) dq + -— ka{q) dq, (36) 



q+^q- """ Jq^^q+ 



where integrals f and f are the limits of the line integrals along 

the segments [g+ — e + ie, q_ + e + ie] and [g- — e — ie, g+ + e — ie] of the 
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contour F/j^, respectively. In view of Lemma [2]^ii) it follows that the last two 
integrals in (15B]1 are equal to 



27rz Jg^ q[q + a) [p2q + /iP2 - A) 



+ ^V^ / , ^ w '^' ^^.q^^z-(q).. + e^-(.i-l)rfg, (37) 

2vr Jg^ q{q + a) (pag + PP2 - A) 

where z~{q) and z"'"(g) are defined in fl3Tl) . Use of the representation e"+*^ = 
e"(cos6 + isinfe) (for a,b G M) completes the calculation of the contour 
integral of 5'(g)/(g + a). D 

Acknowledgements 

F. Avram and M. Pistorius gratefully acknowledge support from the London 
Mathematical Society, grant # 4416. F. Avram and Z. Palmowski acknowl- 
edge support by POLONIUM no 09158SD. Z. Palmowski acknowledges sup- 
port by KBN 1P03A03128 and NWO 613.000.310. M. Pistorius acknowledges 
support by the Nuffield Foundation NUF/NAL/000761/G. 



References 

[1] Asmussen, S. (2000) Ruin probabilities. World Scientific, Singapore. 

[2] Asmussen, S, Avram, F, Usabel, M. (2002) Erlangian approximations 
for finite-horizon ruin probabilities. Astin Bulletin 32, 267-281. 

[3] Ahlfors, L. V. (1979) Complex Analysis. 3rd Edition. McGraw-Hill. 

[4] Bertoin, J. (1996) Levy processes. Cambridge University Press. 

[5] Bertoin, J. (1997) Exponential decay and ergodicity of completely asym- 
metric Levy processes in a finite interval. Ann. Appl. Probab. 7(1), 156- 
169. 

[6] Pervozvansky, A. A. (1998) Equation for survival probability in a finite 
time interval in case of non-zero real interest force . Insurance Math. 
Econom. 23, 287-295. 



18 



[7] Suprun, V.N. (1976) Problem of destruction and resolvent of a termi- 
nating process with independent increments. Ukrainian Math. J. 28, 
39-45. 

[8] Rolski, T., Schmidli, H., Schmidt, V. and Teugels, J. (1999) Stochastic 
processes for insurance and finance. John Wiley and Sons, Inc., New 
York. 

[9] Widder, D.V. (1941) The Laplace Transform. Princeton, NJ: Princeton 
University Press. 



19 



